function [ d ] = rel_compare_profiles( mP, vP, mQ, vQ )
%REL_COMPARE_PROFILES Compute the KL divergence between the two Gaussian
%distributions P and Q

d = 0.5 * (2 * log(sqrt(vQ/vP)) + vP/vQ + (mQ-mP)*(mQ-mP)/vQ -1);

end

